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Abstract 

We examine the effects of pseudoscalar and pseudovector coupling of the tt and r\ mesons in one- 
boson exchange models of the NN interaction using two approaches: time-ordered perturbation 
theory unitarized with the relativistic Lippmann-Schwinger equation, and a reduced Bethe-Salpeter 
equation approach using the Thompson equation. Contact terms in the one-boson exchange ampli- 
tudes in time-ordered perturbation theory lead naturally to the introduction of s-channel nucleonic 
cutoffs for the interaction, which strongly suppresses the far off-shell behavior of the amplitudes 
in both approaches. Differences between the resulting NN predictions of the various models are 
found to be small, and particularly so when coupling constants of the other mesons are readjusted 
within reasonable limits. 

PACS numbers: 13.75. Cs; 21.30.-x; 21.30.Cb; 12.30.Fe 
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I. INTRODUCTION 



Since the discovery and identification of the pion |], as the strongly interacting meson 
anticipated by Yukawa || in 1935, most theoretical efforts to construct a quantitatively 
accurate model of the nucleon-nucleon (NN) interaction have used the pion-nucleon inter- 
action as the first building block. Indeed, phase shift analyses of NN scattering data since 
1959 Q use the one-pion exchange amplitude to fix the phase shifts of the high orbital 
angular momentum partial waves, which are not individually adjusted to fit the data. 

Almost at the beginning of these efforts the question arose of whether the fundamental 
coupling of the pion to the nucleon is of the pseudovector or pseudoscalar type. The question 
arises, of course, from the fact that the fully on-shell one-pion exchange amplitude derived 
from a ttNN interaction Lagrangian with pseudovector coupling is identical in form to that 
from one with pseudoscalar coupling. In early attempts to go beyond the one-pion exchange 
to the two-pion exchange P, pseudoscalar coupling appeared to demand suppression 
of "pair terms" — terms describing the contribution from intermediate states with one or 
more antinucleons. For the exchange of pions with low momenta, this effectively reduced to 
pseudovector coupling, for which the renormalizability of the theory is doubtful, at best |7|, 
||. Pion-nucleon scattering, through the smallness of the scattering lengths, also strongly 
suggested pseudovector coupling. Dispersion-theoretic results for the two-pion exchange 
contribution || [TO]] based on unitarity and analytic continuation of tcN amplitudes to the 
7T7T =f± iViV channel also implicitly favored pseudovector coupling. 

We realize now that the meson theory of nuclear forces is not a fundamental theory 
but, at best, an effective theory. Thus lack of renormalizability in the usual sense is not 
a relevant criterion for the rejection of one form of coupling or another. Furthermore, the 
approximate chiral symmetry exhibited by QCD and implemented with considerable success 
in chiral perturbation theory (%PT) ||11|| , makes it clear that the effective coupling of pions 
to nucleons, at least at low energies, is pseudovector in character. In other phenomena 
the picture is more complex. In pion electro- and photoproduction analyses at low energies 
pseudovector coupling is preferred, whereas at higher energies pseudoscalar coupling provides 
a more economical description [JT2], In the light of this evidence it is clear that any of 
the NN interaction models claiming to be realistic should, to some degree, include pion 
exchange with pseudovector coupling. 

With the exception of some recent models based on baryon x?T [14, 15, 16, 17, lq| , NN 
models of the past three decades include, in addition to the one-pion exchange, contributions 
due to exchange of heavier mesons, whether explicitly, as in the various one-boson exchange 
(OBE) models | _ 



2C, El, 22, 



24, 25, 26, 27, 28 



or implicitly, as resonant t-channel 
exchange of two pions, as in the dispersion-theoretic approaches [§, |l(J. Even within the 
OBE models there is no single preferred approach. There are those models that are based 
on a Bethe-Salpeter [2{| approach and in which the unitarizing equation is some three- 
dimensional reduction of the Bethe-Salpeter equation |22], |24|, |25|, |26|, |30], [31], |32| , and there 
are those, including some of the various Bonn potentials [^, 2£| |3l|, that are based on 



time-ordered perturbation theory (TOPT). Our main interest in this work is in TOPT, but 
we shall also examine differences between pseudoscalar and pseudovector coupling in both 
of these approaches. For that purpose we will utilize a specific three-dimensional reduction 
of the Bethe-Salpeter equation known as the Thompson equation [31]. 

In covariant perturbation theory one starts with a Lorentz-invariant Lagrangian density, 
from which one derives the Hamiltonian density and, from that, the Hamiltonian. For 
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Lagrangians with scalar or pseudoscalar mesons without derivative coupling, the interaction 
part of the Hamiltonian density is just the negative of the interaction part of the Lagrangian 
density. For Lagrangians with derivative coupling or with vector mesons, however, non- 
covariant "contact" terms arise in the Hamiltonian density. These terms are necessary to 
cancel the non-covariant terms in the meson propagators so that, in any order of perturbation 



theory, the resulting amplitude is covariant p£l. From a procedural point of view, this means 



that in the Feynman rules one simply drops the contact terms and the non-covariant parts 
of the propagators. 

In TOPT, however, one does not use particle propagators. They are effectively supplied 
by the vertex functions and energy denominators in the time-ordered diagrams. In order to 
obtain covariant results in TOPT starting from a Lagrangian density with derivative cou- 
pling, or with vector mesons, one must include the contributions of the contact interactions 
in the Hamiltonian in the appropriate order of the perturbation expansion. Therefore, for 
single pion exchange with pseudovector coupling in NN scattering, i.e., in second order in 
the coupling constant, one must include not just the meson exchange diagrams, but also the 
four-point NN contact interaction, as shown in Fig. 1, and similarly for vector meson ex- 
change. Only then will the result agree with covariant perturbation theory when all external 
particles are on their mass shells. 

The main focus of the work that follows is to compare the results of inclusion of the full 
pseudovector-coupled pion exchange with those of pseudoscalar coupling in one-boson ex- 
change models based on time-ordered perturbation theory. This is not simply a moot point: 
the pion-nucleon interaction used in Ref. |33|, although nominally of pseudovector type, is in 
fact pseudoscalar, and thus has a very different off-energy shell behavior from pseudovector 
coupling. Furthermore, contact terms in the vector meson exchange contributions, as well 
as gauge terms, which have heretofore been ignored, will be retained. 

We shall, in this work, restrict ourselves to examining one-boson exchange models, since 
the contact terms in higher orders of TOPT present additional difficulties that are not easily 
resolved. We will, as in all models of this kind, need to introduce cutoffs in order to ensure 
convergence of the integral equations used to unitarize the scattering amplitude. We will 
also compare the results to those obtained from a reduced Bethe-Salpeter approach, in which 
the problem of contact terms doesn't arise. 

In Sect. II we describe the models that we employ in the present study. In particular, we 
give the equations that are used for unitarizing the scattering amplitude and we discuss some 
specific aspects of the vertex form factors. Finally, we outline the strategy that was followed 
for fixing the free parameters of the models. Results for phase shifts as well as for the 
deuteron properties are presented in Sect. III. We analyze the consequences of pseudovector 
versus pseudoscalar coupling in TOPT as well as in a model based on the Thompson equation 
and we also compare the different approaches. In Sect. IV we summarize our results and 
draw some conclusions. Technical details such as the underlying Lagrangians, and the 
potential matrix elements in TOPT and for the Thompson equation are summarized in the 
Appendices A, B, and C. 



II. THE MODELS 

In order to make the comparison of different treatments of pseudoscalar meson exchange 
explicit we will focus on four one-boson exchange models of the NN interaction. The first 
model that we will consider is that of pseudovector coupling of the pseudoscalar mesons, 7r 
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and T], in TOPT. The second will be the same with the exception that pseudoscalar coupling 
of the pseudoscalar mesons will be used. For the third and fourth models we make the same 
comparison of pseudovector and pseudoscalar coupling for the pseudoscalar mesons, but 



within a Bethe-Salpeter approach to the problem using the Thompson equation |31|] to 
unitarize the scattering amplitude. 

The roster of exchanged mesons in each of the four models is identical: pseudoscalar 
mesons 7r and 77, vector mesons p and u, and scalar mesons a and a . Contact terms as 
well as gauge terms arising in the polarization sums in vector meson exchange in TOPT 
will be retained. In the reduced Bethe-Salpeter approach, gauge terms in the vector meson 
propagators will also be retained. The masses and quantum numbers of the mesons are 
given in Table [j| 

The interaction Lagrangian densities for all of the meson-baryon interactions in our model 
calculations are given in Appendix A, along with the corresponding Hamiltonian densities. 
The matrix elements of the corresponding second order potentials for the TOPT-based 
models are presented in Appendix B. Shown schematically in Fig. 2 are the interaction and 
the kinematics for the potential, V, in the NN center-of-mass frame. For this we unitarize 
the scattering amplitude through use of the Lippmann-Schwinger equation, 



V + V T (1) 

W - H n +ie y ' 



or, more precisely, 



< I T(W) I p\iX 2 > = < p'KK I V(W) I p, X1X2 > 

■ < I T(W) I pXrX 2 > . (2) 

Here W represents the total energy of the NN system. 

The Thompson equation, which we use to unitarize the scattering amplitude in the other 
two models, is given by 

< p'X[X' 2 I T(W) I pXiX 2 > = <p'X[X' 2 I V(W) I p, AiA 2 > 

2 

Tfi 

■ < quip* I T{W) I pAiA 2 > . (3) 

The simple prescription for obtaining the matrix elements V from the matrix elements V^W 7 ) 
of Eq. (Q) is given in Appendix C. 

Except for the cutoff functions Fj(W,p',p) that multiply the field-theoretic meson ex- 
change amplitudes in the potential matrix elements V and V, the input to our model calcula- 
tions is completely specified. The cutoff functions, as already mentioned, are needed for the 
convergence of the scattering equation, but the form that one chooses is largely arbitrary. 
One commonly used, especially in one-boson exchange models, is the so-called "multipole" 
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form, 



Fj(W,p',p] 



A 2 - m 2 



J 3 



L A 2 + {p'-p) 2 



(4) 



where Aj and rij are the free parameters and rrij stands for the meson mass. Part of 
the appeal of this particular cutoff is that the connection between the range at which the 
cutoff becomes effective for that meson potential is simply related to the cutoff parameter: 
Rj,cutoS ~ v/tiJ/Aj. Considered as a product of two mesonic form factors, this is then 
interpreted as, or assumed to be, a reflection of the effective size of the meson cloud in that 
part of the NN interaction. 

In treating contact terms in the interaction, however, this form immediately raises some 
difficulties. The first is that, through the dependence on the three-momentum transfer 
in the denominator, it introduces effects of the contact terms into all angular momentum 
partial waves, whereas the contact terms alone, which are polynomials of low degree in the 
sine or cosine of the cm. scattering angle, contribute only to a few of the lowest partial 
wave amplitudes. The contact term arising from pseudovector pion coupling, for example, 
contributes only to s- and p-waves. In adopting the multipole form, one singles out the 
three-momentum transfer as the variable in which to cut off the potential, but there is no 
compelling reason to do so, and good reason not to. Indeed, in the effective field theoretic 
approach of Epelbaum et al. []17|] , a purely s-channel cutoff is used. 

On the other hand, we wish to have the comparison presented in this work make some 
contact with models in the literature [27], that employ t-channel cutoffs. Furthermore, 



we wish to apply the same cutoff to both the meson exchange terms and the contact terms 
in order that one term not receive excessive weight from one kinematic region in the range 
of the loop integral as compared with the other. This is a problem of long standing and we 
will not attempt to address it here. Rather, we will simply adopt ad hoc the form 



Fj(W,p',p[ 



A 2 -m 2 



A? + (P ' - P. 



A 4 \ / A 4 



a% + (w 2 /a-e 2 p ) 2 ) \A% + (wy4-E 2 pl y 



(5) 



that is, we take the form of Eq. (g) and multiply it by a factor A^/(A^ + (W 2 /A - E 2 ) 2 ) for 
each nucleon line with momentum q entering or leaving the interaction. In all cases except 
as noted in Table |T| below we take rij = 2. 

This form provides the potential with both t- and s-channel cutoffs. Although it does 
not remove the objection, raised above, of mixing contact terms into higher partial waves, 
we shall mitigate that effect via the proviso that the parameter Aj for meson terms whose 
interactions include contact terms be chosen large enough so as to have negligible effects in 
partial waves with orbital angular momentum / > 2 , except for 3 Di because of its coupling 
to 3 ^. 

As there are free and nearly-free parameters in the models, a simple comparison of the 
AW phase shifts produced by the four models for a given set of parameters would not, in 
our opinion, be a useful way to present the results of our investigation. Such a comparison 
might be seen to favor one model over another, which is not our intention. 

Instead, for each of the four models, we perform a constrained least-squares fit of the 
adjustable parameters to the AW phase shifts in the range of laboratory kinetic energies 
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from 20 to 300 MeV taken from the energy independent NN phase shift analysis SP40 of 



Ref. 35. By "constrained" we mean that the parameters that are allowed to vary, such 
as coupling constants and cutoff masses, are restricted to a range that is consistent with 
values used in the various meson-exchange models of the NN interaction in the literature. 
The results of these best fits will then be compared in detail. We will also, for the sake of 
completeness, make a comparison of pseudovector and pseudoscalar coupling in the TOPT 
models with identical parameter sets. 

Admittedly, this is not a rigorous procedure, but it reflects better the intention of our 
work. We wish to examine whether the terms that, in principle, should be included in one- 
boson exchange models based on TOPT but have heretofore been omitted require a major 
reworking of previous models, or if their effects can be compensated for by relatively small 
adjustments in the parameters of the other models. Therefore we are not concerned with 
"high precision" fits of the models to the phase shifts, but with qualitatively acceptable fits 
of the magnitudes and energy dependencies of the model results to the data, since further 
refinements of the models, such as the inclusion of two-pion exchange or effects of baryon 
resonances |33|, would necessitate refitting of the parameters and, presumably, result in 
quantitatively better fits. 

At this point we must inject a word of caution for the reader. As our models have 
both s- and t-channel cutoffs to regulate the integral equations, one should not expect the 
cutoff masses that we use, especially those for the mesonic form factors, to agree well with 
those in the literature for models using different cutoff schemes. The way that cutoffs are 
implemented in any model is part of the model and thus have a large influence on the values 
of the model's parameters. 

Our procedure is first to select the parameters in the models that we wish to vary and then 
set the limits of variation of each of these parameter to values that we consider reasonable. 
For example the nN coupling constant, g^/^ir is fixed at 13.8, whereas the rjN coupling 
constant, g^/^ir, since it is less well-known, is allowed to vary from to 6. Cutoff masses, 
Aj, are allowed a fairly large range, but are required to be greater than 1 GeV. The meson 
masses are fixed at the values given in Table |. For the nucleonic cutoff mass, An, we 
explored a range of values and found that good fits with all the models could be obtained 
for Ajy between 600 and 900 MeV. For values below 600 MeV the potential was too strongly 
suppressed and for values above 1 GeV the cutoff had almost no effect. We therefore fixed 
the value of A^r at 700 MeV in all the models. With this choice the contributions to the 
scattering amplitude at low cm. energies from intermediate states with energies above the 
pion production threshold are strongly suppressed. (We should remark here that this value 
of A at is of the same order of magnitude as the s-channel cutoff employed by Epelbaum et 



al. [17] in their effective field-theoretic approach.) We then perform a least squares search 



on the variable parameters. 



III. RESULTS 



The parameters for our four models are shown in Table |J, where we give the complete set 
of coupling constants and cutoff masses for the TOPT models, and in Table pTlj , where we 



show them for the models based on the Thompson equation |3T[|. The phase shifts predicted 
by the various models we have considered are shown in Fig. 3, where we present the results for 
the TOPT models, and in Fig. 4, where we present them for the Thompson equation models 
and, for purposes of comparison, also for the TOPT model with pseudovector coupling. 
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Phase shifts are shown for partial waves only up to J = 3, omitting €3, since the differences 
between the phase shifts calculated with the various models are almost imperceptible in the 
omitted phase shifts. Since the pion coupling is the same in all models, it is clear that the 
higher partial waves in the energy regime considered are practically identical. The deuteron 
properties calculated from the models are given in Table IV. 

As our primary interest is in the differences due to alternative couplings of the pseu- 
doscalar mesons, we shall first examine the results for the TOPT models. We shall then 
turn our attention to the results of Thompson equation models. Afterwards, we shall briefly 
compare the results of the models in the two approaches. 

A. TOPT models 

The two "best fit" TOPT models — curves PV and PS in Fig. 3 — both give reasonably 
good descriptions of the phase shift data. Where there are discernible differences, it is 
difficult to form a consistent picture of the effects of the two alternative couplings that 
cannot be compensated by relatively small adjustments in the cutoffs or in the coupling 
constants of the other mesons. As mentioned before, the pion coupling constant is the same 
in all models. In order to demonstrate the difference of the two coupling schemes as it arises 
when there is no readjustment of the other meson parameters, Fig. 3 includes a calculation, 
labelled PSO, in which the only difference from model PV is the change of the coupling of the 
71 and i] mesons from pseudovector to pseudoscalar. A comparison of the results of PV with 
PSO shows that, even with no readjustments of the cutoffs or coupling constants of the other 
mesons, the differences between pseudovector and pseudoscalar coupling for n and rj is quite 
small, especially at low energy. There is, however, one glaring exception: the predictions 
for the mixing parameter e\. This parameter as calculated from PSO is almost double the 
value from PV — and the data — throughout the energy range shown, thus indicating that 
the tensor force resulting from the two coupling schemes is quite different. 

Indeed, it is the need to describe e\ more accurately that largely drives the changes in 
the coupling constants and cutoffs in model PS from model PV. With very few exceptions, 
the results of PSO are closer than those of PS to the results of PV. The exceptions, as one 
might expect from their connection with e±, are the phase shifts 3 Si and 3 Di, to which e\ is 
coupled. 

Comparing the coupling constants and cutoffs of PV and PS, one sees that the r\ is 
essentially unchanged. The p coupling in PS is slightly smaller, but that is partially offset 
by a cutoff mass that is slightly larger. The greatest differences between the parameters 
in the two models are in the a and uj. The u coupling in PS is about 25 percent larger 
than in PV, although its cutoff mass is smaller, which tends to compensate for the increased 
repulsive strength at short distance. The increase in the a coupling is necessary, apparently, 
to provide attraction at intermediate range to counter the greater repulsion due to the u. 

While the picture isn't entirely clear, the competition between increased attraction due 
to the larger a contribution at intermediate range and increased repulsion at short range 
due to the larger u contribution appears evident. In the mid-peripheral uncoupled phases 
where differences can be observed, i.e. and 3 D2, PS is slightly more repulsive than PV, 
while in the more peripheral phase shift ^3, the reverse is true. For higher partial waves, 
which are not shown in the figures, the phases are almost entirely given by the on-shell pion 
exchange term, but the effect is, as one would expect, the same as in 1 F 3 , although the 
differences are so small as to be invisible on graphs. 
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B. Thompson equation models 



The phase shift fits for the two models that employ the Thompson equation are shown 
in Fig. 4, along with the results for the TOPT-based model PV discussed in the previous 
section. The curves labelled TPV show the results for the model with pseudovector cou- 
pling of the 7r and r\ mesons, while those labelled TPS show the results for the model with 
pseudoscalar coupling of the 7r and rj mesons. Both TPV and TPS are best-fit results in 
the restricted sense discussed in Section II, in which bounds are placed on the range of the 
adjustable parameters. 

As in the case of the TOPT models, the phase parameter that shows the most striking 
difference between pseudovector and pseudoscalar coupling is the mixing parameter e±. Using 
the models with the Thompson equation with the cutoffs implemented as we have described, 
it is impossible to achieve a satisfactory description of ei while simultaneously keeping the 
coupling constants within reasonable bounds when pseudoscalar coupling is used for the 7r 
and rj mesons. 

Apart from ei, the only other phase shifts that reveal any noticeable difference between 
pseudovector and pseudoscalar coupling in the Thompson equations models are 1 Pi, 62, and 
3 D3, and even there the differences are rather small. For the most part, the two Thompson 
equation results are closer to each other than either is to PV or to PS. This result isn't 
surprising, since most of the parameters in TPV and TPS are the same, which reflects the 
fact that they are at the limits of their permitted ranges. An interesting result is that the 
t-channel cutoff of the a meson, A CT , is extremely large — 10 GeV, which is effectively infinite. 
This suggests that the s-channel cutoff has a very powerful effect in the Thompson equation 
models, and the large value of A CT reflects an effort of the fitting program to increase the 
attractive effect of the a exchange contribution to counter the increased strength of the u 
contribution. Indeed, with very few exceptions, both of the Thompson models are more 
repulsive than PV, as one might expect from the large value of g%. 

A slightly different view of the models that we have considered is provided by the deuteron 
parameters compiled in Table ffVL which shows some small — but consistent — differences be- 
tween them. All four models are adjusted to fit the deuteron binding energy very accurately. 
Both TOPT models give a very good value of the quadrupole moment with a relatively low 
d-state probability, which is characteristic of TOPT models, although the tendency of the 
Thompson models to have relatively large d-state probability is mitigated here, presumably 
because of the strong s-channel cutoff that was not present in earlier work [T2(J. In both 
the TOPT and Thompson equation models, pseudoscalar coupling results in a lower d-state 
probability than pseudovector coupling. The asymptotic s-wave, As, is somewhat high for 
the TOPT models, although the asymptotic d-to-s ratio, Ap/As, is acceptable. As expected, 
the d-state probability of the Thompson models is larger than in the corresponding TOPT 
models, resulting in a lower value of Ag. The value of Ajj/Ag, however, is not very different 
among the four models, which probably explains the smaller quadrupole moments for the 
Thompson equation models: the slightly greater d-state probability in the Thompson equa- 
tion deuteron wave functions is not sufficient to counter the slightly more compact structure 
of the deuteron that the Thompson equation produces. 
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IV. SUMMARY AND CONCLUSIONS 



It is clearly impossible to draw any general conclusions concerning the effects of pseu- 
doscalar as opposed to pseudovector coupling of pseudoscalar mesons in the NN system. 
Any statements that we make are necessarily qualified by their model dependence. It is 
nevertheless useful to summarize our approach and findings and to note, if possible, any 
tendencies within the limited context of the models that we have investigated. 

In the first place, we wished to treat pseudovector coupling properly in the TOPT ap- 
proach, using the Lippmann-Schwinger equation to generate unitary scattering amplitudes 
from the one-meson exchange amplitudes, and to compare it with pseudoscalar coupling. 
The presence of contact terms in the pseudoscalar meson exchange terms, as well as in the 
vector meson exchange terms with tensor coupling, led us to introduce nucleonic — i.e. s- 
channel — form factors in addition to the t-channel form factors that are usually employed 
in meson exchange modes of the NN interaction. 

For the purpose of comparison with a different approach, we examined the difference 
between the two coupling schemes in the context of a model based on a particular version 
of three-dimensional reduction of the Bethe-Salpeter equation, the Thompson equation. In 
this approach the energy denominators and the off-shell continuations of the meson exchange 
amplitudes differ from those of TOPT. In particular, no contact terms appear. In order to 
keep the comparison as close as possible, we chose similar form factors to those in the 
TOPT-based models that we studied. 

Within each model we made a restricted best fit to the NN data, allowing meson coupling 
constants and form factors to vary within broad limits chosen with regard to values of these 
parameters found in earlier works. We also examined the effect of simply changing the 
pseudovector coupling of the n and 77 to pseudoscalar, within the TOPT approach. 

Not surprisingly, perhaps, most of the differences in the scattering and bound state prop- 
erties calculated with the various models were quite small, which probably reflects the fact 
that the cutoffs were strong enough to strongly suppress differences in the off-shell behavior 
of the two coupling schemes. The most dramatic difference in the TOPT models appeared 
in the mixing parameter e\ when the simple change of pseudovector to pseudoscalar coupling 
was made. This should not be surprising, since the effect of the contact terms is limited to 
s- and p-waves and states coupled to them. The diagonal pion exchange amplitude is rather 
weak so that it is in the relatively small mixing parameter e±, which is dominated by 7r and 
p exchange, that changes in the short-range behavior of the tensor force is most strongly 
felt. Indeed, most of the differences between the "best fit" pseudovector and pseudoscalar 
models was due to the readjustment of parameters in the latter needed to produce a better 
fit to e\. 

The same effect was observed in the Thompson equation models, with the largest rela- 
tive difference between the pseudovector and pseudoscalar coupling versions remaining after 
refitting coupling constants appearing in e\. Other differences between the two Thompson 
equation models were very small, owing in large part to the fact that several of the adjustable 
parameters were at the limits of their permitted ranges in both cases. 

Differences in the deuteron parameters among the models were similarly quite small, with 
only the slight tendencies in d-state probability, quadrupole moment and asymptotic s-state 
predictions noted previously. 

The motivation of this study was to see whether, with a very restricted set of models, one 
type of coupling of the pseudoscalar mesons would be able better to reproduce NN scatter- 
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ing data than the other. Within the TOPT approach one might claim that, on the whole, 
pseudovector coupling yields a slightly better description of the data than pseudoscalar, but 
not in every partial wave. Within the Thompson models, the differences between pseudovec- 
tor and pseudoscalar coupling are still smaller. Between the TOPT and Thompson model 
predictions there is no clear best. We must conclude, therefore, that the results obtained 
for iViV scattering with the models considered in this work present no compelling evidence 
that one form of coupling of the pseudoscalar mesons to nucleons is favored over the other. 
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APPENDIX A: INTERACTION LAGRANGIANS AND HAMILTON! ANS 



We present here, for the purpose of completeness, the interaction Lagrangian densities 
and the corresponding interaction Hamiltonian densities for the various meson-nucleon in- 
teractions used in this work. 



1. Scalar meson 



£i, s = -9s^4>s (Al) 
n I>s = -d, a (A2) 



2. Pseudoscalar meson, pseudoscalar coupling 



£i, PP s = -ig P ipl 5 ip(/>p (A3) 
Hi,pps £*i,pps (A4) 



3. Pseudoscalar meson, pseudovector coupling 



Ci, PP v = -^^V^0 P (A5) 
Hj, ppv = -L lwo + i^^yVV (A6) 



4. Vector meson 

Here we take for the free meson Lagrangian the form 



d,v = -\f»„F^ + l -mlA^, (A7) 
where = d )X A v — d u A^. We then have 

Ci, v = -grd-f^Ap - JjL-ja^F^ (A8) 
^ = -^ + ^(^) 2 + I^(^ )2 (A9) 

where a^ v = |[7 M , T^]-- We use the Bjorken-Drell conventions, summing on repeated indices 
and using the Latin letter % to denote spatial indices 1..3. Isotopic spin notation has been 
suppressed. 
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APPENDIX B: MATRIX ELEMENTS OF THE POTENTIAL 



We show here the potential matrix elements in time-ordered perturbation theory to second 
order in the meson-nucleon coupling constants derived from the Hamiltonians in Appendix 
A. The field-theoretic matrix element for the contribution of meson j is multiplied by a 
cutoff function Fj(W,p',p), as described in Section II. 

We use the helicity basis in the two-nucleon center-of-mass (cm.) frame. The total cm 
energy of the system is W, the momenta of the ingoing and outgoing nucleons are (p, —p) 
and (p 1 , —p') with corresponding energies E p and E p > and helicities Ai, A2 and A^, A2, 
with A = ±±. The energy transfer in the interaction is S = E v < — E p , the three-momentum 

~ — * 

transfer in the interaction is k = p — p ' and the energy of the exchanged meson of type 
j is uj J k = \Jk 2 + m 2 -. The energy denominator for the meson exchange term is given by 
D> = W-E P - E p , - u{. 

1. Scalar meson 



k (2) = ^^{[iiitiwwr.Pi (bi) 



2. Pseudoscalar meson, pseudoscalar coupling 



= -jif l ^{bWh} F P(PS) (W,P',P) (B2) 



3. Pseudoscalar meson, pseudovector coupling 



(, ,P\2 _ A2 \ 1 



+ Li [r7°]i[7 5 7°]2 ) - -^[7 5 7°]i[7 5 7°] 2 \F p(j)v) (W,p> ,p) (B3) 
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4. Vector meson 



+ -pa)/j]2 + K\p\-Pi) a W%) + (K) 2 - S 2 )[a%[a°% 

+ (2^ 3 2m^ Il2 ^{^ ([7ll[ ^ (P2 " P2)ah + P*(A-Pi)a]i[f\2) 
+ 5([YW% + MiM^.W^,^ (B4) 
In the expressions above we have used the compressed notation 

Vf } =< p'X[X 2 I vj 2) (W) I p\i\ 2 > (B5) 

and 

[A] 1 [B] 2 = [u(p', X[)Au(p, X,)] [u(-p>, X' 2 )Bu(-p, A 2 )] (B6) 

for the matrix elements of the potential. 

The isospin factor I12 is 1 or f\ ■ f 2 as the isospin of the meson concerned is or 1. The 
normalization of the Dirac spinors is 

^(f,X)u(p,X) = l, (B7) 

and the four- vectors p\ and p 2 are (E p ,p) and (E p , —p), respectively, and similarly for the 
primed quantities. 

We have written the matrix elements in a way to distinguish clearly the meson exchange 
terms, which contain energy denominators = W — E p > — E p — ou J k , from the contact terms, 
which do not. The equivalence of the fully on-shell one-pion exchange with pseudovector 
coupling with that of the one with pseudoscalar coupling is then evident. In that case 
E p = E p i = W/2, D p = —ijj v k and 5 = 0. The term containing {oJ p k ) 2 is exactly canceled by 
the contact term. 

APPENDIX C: THE THOMPSON EQUATION 

Matrix elements for the potential to be used in the Thompson equation can easily be 
found using the results of Appendix B by applying the following prescription: 
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1. Change the nucleoli spinor normalization to 

u(p,\)u(p,\) = l. (CI) 



2. For each meson contribution Vj 2 \ replace the energy denominator w^Di accord- 
ing to 

ujpi — -K) 2 = -(?' - Pf - ™] (C2) 



3. Drop all contact terms. 
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TABLE I: Quantum numbers and masses of the mesons used in the models in this work. 



meson I (J p ) mass (GeV) 

7T 1 (0~) 0.13803 

rj (0~) 0.5488 

cr (0+) 0.52 

a 1 (0+) 0.983 

p 1 (1~) 0.769 

u 0(0") 0.7826 



TABLE II: Parameters for the TOPT models. The columns labelled PV and PS represent, 
respectively, the values for the model with pseudovector and pseudoscalar coupling of the n and 
rj mesons, as in the figures. Fixed parameters are shown in parentheses. The values of searched 



parameters are determined by a least-squares fit to the SP40 AW phase shift analysis of Ref. 35. 



For these models n n 


= 1 (see Eq. (|)). 










Coupling constant 


PV 


PS 


Cutoff 


PV(GeV) 


PS(GeV) 




(13.8) 


(13.8) 


K 


2.50 


1.80 




2.15 


2.00 


A,, 


1.00 


1.21 


all^ 


6.443 


7.177 


A. 


10.00 


9.80 


alj^ 


0.858 


1.859 




2.50 


1.50 




1.14 


1.06 


Ap 


1.45 


1.73 


fp/dp 


5.12 


4.40 








all^ 


17.40 


22.20 


A w 


1.41 


1.35 


fu/gu 


(0) 


(0) 














A N a 


(0.7) 


(0.7) 



a See Section II, Eq. (5) for the use of A 



AT- 



TABLE III: Parameters for the Thompson equation models. The columns labelled TPV and TPS 
represent, respectively, the values for the model with pseudovector and pseudoscalar coupling of 
the 7r and rj mesons, as in the figures. Fixed parameters are shown in parentheses. The values of 
searched parameters are determined by a least-squares fit to the SP40 AW phase shift analysis of 



Ref. 35 



Coupling constant 


TPV 


TPS 


Cutoff 


TPV(GeV) 


TPS(GeV) 


all^ 


(13.8) 


(13.8) 


A^ 


2.00 


1.80 


aV^ 


4.22 


5.00 


A^ 


1.00 


1.00 


6.871 


7.205 


K 


10.00 


10.00 


alj^ 


5.00 


5.00 


Aa 


2.50 


1.50 




0.800 


0.800 


Ap 


1.34 


1.31 


fp/ap 


6.89 


6.89 










25.00 


25.00 




1.245 


1.260 


fcj/av 


(0) 


(0) 









A N (0.7) (0.7) 
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TABLE IV: Deuteron properties calculated with the four models considered. 



Quantity 


Experiment 


PV 


PS 


TPV 


TPS 


-E d [MeV] 


2.24575(9) |§ 


2.22450 


2.22447 


2.22453 


2.22459 


Pd [%] 




3.8 


3.4 


4.4 


3.9 


Qd [fm 2 ]° 


0.2859(3) |7|, 38 


0.2784 


0.2780 


0.2765 


0.2754 


A s [fm- 1 / 2 ] 


0.8846(9) |7|, |39 


0.9117 


0.9119 


0.8974 


0.8950 


A D /A s 


0.0256(4) || 


0.0255 


0.0260 


0.0252 


0.0257 



"Theoretical values do not include meson exchange current contributions. 
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FIG. 1: Schematically, the three terms in time-ordered perturbation theory that may contribute 
to NN scattering in second order of the coupling constant in a meson exchange model. 




P' E P -P' E P 



FIG. 2: . Kinematics for the NN potential in the NN center- of-mass frame with total energy W. 
The fermion lines are labelled by their 3-momentum, energy and helicity. 
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FIG. 3: . Phase shifts for TOPT models. The curves labelled PV are the results of the best fit 
for the model with pseudovector coupling of the ir and the 77; those labelled PS are the best fit for 
pseudoscalar coupling. The curves labelled PSO are the results of using the parameter set of PV 
but with pseudoscalar coupling. Parameters for the models are given in Table |fij. The triangles 
represent the phase shift analysis by the Nijmegen Group at selected energies Q, and the open 
circles stand for the energy independent analysis SP40 from the CNS DAC Services SAID |35|| . 



19 





20 




Fig. 3 cont'd 
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FIG. 4: . Phase shifts for Thompson equation models. The curves labelled TPV are the results of 
the best fit for the model with pseudovector coupling of the tt and the 77; those labelled TPS are 
the best fit for pseudoscalar coupling. Included for purposes of comparison are the results of the 
best fit of the TOPT model with pseudovector coupling (curves labelled PV). Parameters for the 
models are given in Table |D|. The error bars are the same as in Fig. 3 
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Fig. 4 cont'd 
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